Instantons, diquarks and non-leptonic weak decays of hyperons 



o 
o 

(N 

wo: 
< 

(N 
> 
O 
00 

(N 
O 
^t" 
O 

^ 

a: 

CD 
43 



X 



M. Cristoforetti 1 , P. Faccioli 23 , E.V. Shuryak 4 and M. Traini 35 

Dipartimento di Fisica, Universita di Milano 
2 ECT* , Villazzano (Trento) 
3 INFN, Gruppo Collegato di Trento 
4 Department of Physics and Astronomy, 
SUNY at Stony Brook 
5 Dipartimento di Fisica, Universita di Trento. 
(Dated: February 1, 2008) 

This work is devoted to the study of the non-perturbative contributions in non-leptonic hyperon 
decays. We show that the instanton-induced 't Hooft interaction can naturally explain the A I — 1/2 
rule, by generating quark-diquark clustering inside octet baryons. We compute P-wave and S-wave 
amplitudes in the Instanton Liquid Model, and find good agreement with experiment. We propose 
a model-independent procedure to test on the lattice if the leading quark-quark attraction in the 
+ anti-triplet channel responsible for diquark structures in hadrons is originated by the interaction 
generated by quasi-classical fields or it is predominantly due to other perturbative and/or confining 
forces. 



I. INTRODUCTION 

Weak decays of hadrons encode important informa- 
tion about the meson and baryon structure and about 
the QCD interactions in the perturbative and non- 
perturbative regimes. The natural scale of weak pro- 
cesses -set by W boson mass- is much larger than all 
other scales involved in the hadron internal dynamics. 
This implies that weak interactions are effectively local 
and therefore can resolve short distance structures inside 
hadrons. Moreover, their explicit dependence on quark 
flavor and chirality can be exploited to probe the Dirac 
and flavor structure of the non-perturbative QCD inter- 
action. 

Among the large variety of weak hadronic processes, 
a prominent role is played by the non-leptonic decays 
of kaons and hyperons, which are characterized by the 
famous A I — 1/2 rule Q. With this name, one refers 
to the empirical observation that amplitudes in which 
the total isospin is changed by 1/2 units are roughly 20 
times larger than the corresponding amplitudes in which 
the isospin is changed by 3/2 units. 

Despite nearly 40 years of efforts, the microscopic dy- 
namical mechanism responsible for such a striking phe- 
nomenon is still elusive. Neither electro- weak nor per- 
turbative QCD interactions can account for the dramatic 
relative enhancement of the A / = 1/2 decay channels. 
Its origin must therefore reside in the non-perturbative 
sector of QCD. 

Important insight on the role of non-perturbative dy- 
namics in non-leptonic hyperon decays has come from the 
observation that in the pole-model (see below) the sup- 
pression of the decays in the A / = 3/2 channel can be 
explained if the quarks participating to the weak decay 
are in an anti-symmetric color combination (Pati-Woo 
theorem, 0). Unfortunately, in a simple Constituent 
Quark Model picture it is not easy to obtain satisfactory 
quantitative predictions for both the P-wave and the S- 
wave amplitudes. One usually needs to make additional 



model-assumptions on the pole-model part of the am- 
plitude, and this somewhat spoils the simplicity of the 
approach. For example, in order to reproduce the data 
on S-wave amplitudes, one needs to include l/2~ inter- 
mediate states . 

From these considerations it follows that further in- 
vestigations are still needed in order to understand 
the non-perturbative QCD dynamics underlying non- 
leptonic weak decays. In particular, it would be desirable 
to set-up a field theoretic calculation which accounts ex- 
plicitly for the current quark and gluon degrees of free- 
dom. In this work we explore the possibility that the 
phenomenology of hyperon decays can be understood in 
the Instanton Liquid Model (ILM). Such an approach is 
derived directly from the QCD Lagrangian, by selecting 
a specific set of gauge configurations which are assumed 
dominate the path integral. 

Instantons are topological gauge configurations which 
dominate the QCD path integral in the semi-classical 
limit. They generate an effective quark-quark interac- 
tion ('t Hooft vertex) which breaks spontaneously chiral 
symmetry and solves the U(l) problem Evidence for 
instanton-induced dynamics has been accumulated over 
the years from a variety of phenomenological studies [f| 
as well as from lattice simulations [a, l3> 13 13 ■ I n general, 
these non-perturbative vacuum fields play an important 
role in the chiral dynamics of light quarks [Tol ] , but it is 
generally believed that they do not provide an areal law 
for the Wilson loop. 

The ILM assumes that the QCD vacuum is saturated 
by an ensemble of instantons and anti-instantons. The 
two phenomenological parameters of the model are the 
instanton average size (p — 1/3 fm) and average density 



(n 



1 fm ). These values were first extracted 20 



years ago from the global properties of the QCD vacuum 
(quark and gluon condensates) [l2^ . 

In the ILM, quarks are bound by the 't Hooft interac- 
tion. Even in the absence of confinement, the structure 
of the lowest-lying part of the light meson and baryon 
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spectra is very well reproduced 0, Q, 0] . In particu- 
lar, in this model the lightest octet of pseudo-scalar and 
vector mesons, as well as the lightest octet and decu- 
plet of baryons, have very realistic masses. Moreover, 
the short-range forces generated by instantons allow to 
reproduce the available experimental data on the pion 
and nucleon form factors and more generally explain the 
delay of the onset of the asymptotic perturbative regime, 
in hard exclusive reactions 0, 0] . 

Besides providing a successful overall description of the 
light hadron phenomenology, instantons have a specific 
property which makes them natural candidates for the 
solution of the A 7 =1/2 problem. In fact, Stech, Neu- 
bert and Xu pointed out that the body of data on non- 
lcptonic kaon and hyperon decays can be simultaneously 
reproduced, if one assumes that the non-perturbative 
quark-quark interaction in the color anti-triplet chan- 
nel is sufficiently attractive to form colored quasi-bound 
structures (diquarks) inside hadrons |23j . Instantons pro- 
vide a microscopic mechanism which generates such a 
strong attraction binding scalar diquarks |3| and lead- 
ing to quark-diquark clustering inside the octet baryons. 

In the past there have been few attempts to under- 
stand the A 7 =1/2 rule with instantons EUcJ. In 
[l8| Kochelev and Vento (KV) computed the instanton 
contribution to non-leptonic kaon decays. On a quali- 
tative level, they found that the inclusion of the instan- 
ton effects indeed produces a strong enhancement of the 
A I = 1/2 decay channel. On a quantitative level, such 
an enhancement was found to be still insufficient to repro- 
duce the experimental data. However, it should be men- 
tioned that non-leptonic kaon decays in the A 7 = 1/2 
channel receive large contribution also from final-state 
interactions, which have not been included in the KV 
analysis. Moreover, it is now clear that the KV calcula- 
tion is undershooting the instanton contribution [3^. 

In |2(j the instanton-induced corrections to the effec- 
tive Hamiltonian for AS — 1 transitions were analyzed 
in the framework of the Operator Product Expansion 
(OPE). They found that such "hard" instanton effects 
are rather small. This result is not surprising: the in- 
stanton field cannot transfer momenta much larger than 
its inverse size 1/p ~ 0.6 GeV, so instanton effects above 
such a scale are exponentially suppressed. For this rea- 
son, in order to draw conclusions about the role played 
by the 't Hooft interaction in weak decays, one neces- 
sarily needs to include their contribution to the "soft" 
hadronic matrix elements. In view of these arguments, 
in the present analysis we shall neglect all instanton cor- 
rections above the hadronic scale set by the inverse in- 
stanton size p — 1/p and compute their contributions to 
low-energy matrix elements. 

The paper is organized as follows. In section [H] we 
analyze the structure of the effective Hamiltonian for 
A S = 1 transitions and explain in detail why instan- 
tons are expected to produce strong enhancement of the 
matrix elements associated to A 7 = 1/2 transitions. In 
section UlTl we review the framework which allows to con- 



nect parity-conserving and parity-violating decay ampli- 
tudes to low-energy matrix elements of local operators. 
The calculation of the decay amplitudes in the ILM is 
presented in section IIVI In section we discuss our 
results and address the question of how to check our 
model assumption of instanton domination for the light 
hadron dynamics. We shall propose a systematic pro- 
cedure to determine on the lattice if the strong quark- 
quark attractive interaction in the anti-triplet + chan- 
nel ( which drives the A 7 = 1/2 rule ) is predominantly 
due to quasi-classical gauge configurations or is instead 
generated by other non quasi-classical fields, associated 
to quark confinement. All results and conclusions are 
summarized in section IVT1 

II. A S = 1 EFFECTIVE HAMILTONIAN AND 
THE ORIGIN OF THE A I = \ RULE 

To lowest order in the Weinberg-Glashow-Salam La- 
grangian, non-leptonic weak decays are driven by a sin- 
gle W-boson exchange. However, such processes receive 
also QCD and QED corrections. These contributions 
are usually included in the framework of OPE, in which 
one separates short-distance "hard" dynamics from large- 
distance "soft" dynamics. The former interactions can be 
treated perturbatively and give rise to the well-known ef- 
fective weak Hamiltonian, which for A S = 1 transitions 
reads 22]: 

H ef S f 1 = VudVus { E c i ( f x)Q i +h. C .\. (1) 

V i=±,3,5,6 

Gf is the Fermi's constant, V u d and V us are quark mixing 
matrix elements, Qi are local four-quark operators and 
Ci(n) are the corresponding Wilson coefficients (p is the 
hadronic scale). The local operators Qi can be written 
as: 

Q± = ~ [(us) V - A (du) v ^ A ± (ds) V - A (uu) v _ A ] 

Q3,5 = (ds)v-A 2J (qq)vTA 

q—u,d,s 

Qe = -2 (qs)s+p(dq)s-P, (2) 

q— u,d,s 

where we have adopted the notation (q q)v±A = ?7p (1± 
75) Q, and (q q) s ±p = q (1 ± 7s)<7 

For a typical hadronic scale, /i ~ 1 GeV, the numerical 
values of the Wilson coefficients are c+ = 0.72, c_ = 
1.97, c 3 = -0.005, c 5 = 0.003, c 6 = -0.008 From 
these numbers it follows that non-leptonic weak decays 
are driven by the terms proportional to the operators Q + 
and Q_, while all other terms can be neglected. 

It is straightforward to verify that the operator Q_ 
triggers decays with A 7 = 1/2, while the operator Q + 
induces transitions both in the A 7 = 1/2 and in the 
A 7 = 3/2 channel. Hence, in order to explain the 
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A I = 1/2 rule, one needs to understand the dynamical 
mechanism which enhances the contribution of the term 
proportional to Q-. 

Accounting only for weak interactions one finds c + = 
c_ = 1 and C3 = C5 = C6 = 0. Clearly, perturba- 
tive strong forces do indeed provide a relatively small 
enhancement of A I = 1/2 transitions. On the other 
hand, a factor 10 is still missing in order to reproduce 
the experimental data. This must necessarily come from 
the non-perturbative sector of QCD. In the OPE for- 
malism, large-distance strong dynamics enters through 
the low-energy matrix elements of the effective Hamil- 
tonian (JJ. Hence, we conclude that non-leptonic weak 
decays are driven by non-perturbative forces which en- 
hance by roughly one order of magnitude the hadronic 
matrix elements of Q-, relative to the matrix elements 
of Q+. 

Significant progress in trying to understand these non- 
perturbative effects has been made in a series of works 
by Stech, Neubert, Xu and Dosch (SNXD) 0, 0, IU 
|2a,|23|. Their starting point was the observation that the 
effective Hamiltonian could be Fierz-transformed into: 



d V us |c_ (ud)l« {su) 3 * 
+c+(/i)(ud)|(sit) 6 + ... + h.c.}, (3) 



where (su)y = eikj sj C '(1 — 75)1^- is a scalar and pseudo- 
scalar color anti-triplet diquark current, while (su)e is 
the corresponding color sextet current (the other cur- 
rents are given by similar expressions). From © it fol- 
lows immediately that the matrix elements proportional 
to c_ (/i) will be greatly enhanced if the non-perturbative 
quark-quark interaction is very attractive in the color 
anti-triplet channel. This is most evident in hyperons: if 
the non-perturbative forces are so strong to allow -say- a 
s and u valence quarks in a E + to form a + anti-triplet 
quasi-bound state, then these quarks will have a much 
larger chance to be caught in the same point and annihi- 
lated by the local {su)y operator in the effective Hamil- 
tonian. A similar ar gum ent can be formulated also in the 
case of kaon decays |37j . Based on this simple dynamical 
assumption, SNXD proposed a phenomenological model 
which simultaneously explains kaon and hyperon non- 
leptonic decays. 

In order to justify the phenomenological assumptions 
of the SNXD model and make contact with QCD, we need 
to identify some non-perturbative gauge configurations 
which, on the one hand, play an important role in the 
hadron internal dynamics and, on the other hand, gen- 
erate color anti-triplet quasi-bound diquarks. Instantons 
have precisely this property. In 14] it was shown that 
the 't Hooft interaction does indeed form a bound anti- 
triplet scalar diquark of mass of roughly 400 MeV. It is 
therefore natural to ask whether these fields can provide 
the microscopic mechanism underlying the A / = 1/2 
rule. 



FIG. 1: Factorization contribution to the non-leptonic hy- 
peron decay (nB'\H e f f\B) . 



III. LOW-ENERGY MATRIX ELEMENTS 

Non-leptonic hyperon decays can be parametrized in 
terms of two constants corresponding to parity-violating 
and parity-conserving transitions: 



{B'Tv\H eff \B)=iu B ' [-A-B75] u B , 



(4) 



where B (B') denotes the initial (final) baryon, and A 
and B are respectively called S-wave and P-wave ampli- 
tudes. The calculation of these amplitudes is generally 
performed by analyzing separately two different contribu- 
tions which correspond to different mechanisms through 
which the pion in the final state in Q can be produced. 

In the so-called "factorization" part of the ampli- 
tude the final meson is excited directly by the 
color singlet axial-vector current present in the effective 
Hamiltonian (as pictured in Fig. The correspond- 
ing parity-conserving and parity-violating amplitudes for 
non-leptonic hyperon decays with tt~ in the final state 
are l23l: 



.4 



7T (fact) 



(dOO + 2c 6 ( M )-^) F w (Mi - Mj) Fjf 



B 



it (fact) _ 



= - Ci (/x)-2c 6 (m)— )F n (M i +M j 



K 



(5) 



with 



m 



K 



m u + md m s + m u 



m 



K 



f/g I Ivy, 



F v = 132 MeV. (6) 



Decay amplitudes with 7T° in the final state are obtained 
from the substitution 

A jJ (fact) = _ 1 A fact( Cl ^ p4+i5 ^ 

V2 



B 



tt (/act) 



1 



V2 



B fa _ ct { Cl -» -c 2 , F 4+i5 -> F 6+l7 ). 



(7) 



In JSJ) the i and j indices select the baryons in the initial 
and final state, and the constants Fji and Gji are the 
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FIG. 2: Pole contributions to the non-leptonic hyperon decay 
(nB'\H eJf \B). 



axial-vector and vector form factors at zero momentum 
transfer, denned as: 

{B j {l/2+)\j;\B i {l/2+)) k ^ Q = *£u(j)7p«(<) 
(B j (l/2+)\J^\B i (l/2+)) k ^ = G^u^Wuii). 



(8) 



Assuming SU/(3) flavor symmetry and using the 
Goldberger-Treiman relation we have: 

g% = V / 2(* fjai F + dj a i D)g with 
V2 

Fji — i fjai, (9) 



9Ji = -^G^iiMi + Mi) 



Notice that, in the flavor symmetric limit, the factoriza- 
tion part of the amplitudes is completely determined in 
terms of experimentally measured low-energy constants. 
In this work, we use the values 1251: 



g = 13.5, 



F + D = l. 



D 



1.6 



(10) 



where the D/F ratio is extracted from semi-leptonic de- 
cays [28|. 

It is immediate to verify that factorization amplitudes 
alone cannot explain the non-leptonic low-energy decays 
of kaons and hyperons [3^| . 

The leading contribution to such reactions emerges 
from a soft-pion analysis of the matrix element By 
applying the PCAC relation, the pion in the final state 
is replaced by an additional operator, expressing the di- 
vergence of the axial- vector current: 



(B j 7r a (q)\H eff (0)\B i ) = toi 



V2(-q 2 + ml) 



F„ ml 

x / <fixe i **{B j \T(d»JgJx)H eff {0))\Bi). 



(11) 



One then applies the well-known identity 

i f #xT{d*J$^x)H afl {0))e i *> 



g" / d*xT(J§Jx)H eff {0))-i[I$,H eff ], 



(12) 



( If is the axial charge operator ) and performs the an- 
alytic continuation to q^ — ► (soft-pion hypothesis). 

The first term in the right-hand side of Ijl2(l leads to the 
so-called "pole contribution" . Physically, it corresponds 
to the processes in which the effective Hamiltonian mixes 
the initial or final baryon with some intermediate virtual 
state (see Fig. I2J - The final results for the pole contribu- 
tions read: 



(pole) 
j'i 

(pole) 
ji 



V2(Mj + Mi 



F„ 



h n Gu 



(Mi - Mi) (Mj - M t ) 



V2 



E 3l h U ~ h jl E H 



(13) 



where M* denotes the masses of the intermediate 
Bi(l/2~) baryon which is mixed with the l/2 + baryon 
by the effective Hamiltonian. The low-energy constants 

hti and E% are defined as: 

J l> J L 

{B j (l/2+)\H% f \B i (l/2+)) = h±u(3)u(i) 
{B j {l/2+)\H% f \B i {l/2-)) = hJ^UMt), 



(B j (l/2-MJB i (l/2+)) k ^ = £S5(i)7M«(i). 



(14) 



i e ff is the parity-conserving part of the effective 
Hamiltonian, and reads: 



where 



Hl° ff = A [eijk(diCuj)ei m h(diCu m ) + 
e ljk {d l C j 5 Uj)eimk(diC 7 5 u m )], 

G F 

A = —= sin 9 C cos 9 cC-(ijl). 
v2 



(15) 



(16) 



H^ff is the parity-violating part of the effective Hamil- 
tonian and reads: 



TJPV _ 



-A[eijk(diCuj)ei mk (diC j 5 u m ) + 
eijk(d~iC j 5 Uj)ei mk (diC u m )], 



(17) 



Using SU/(3) symmetry, one can express these matrix 
elements in terms of few coefficients: 
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hj = eS j6 

Eft = 2V2 (if jai F- 

Eni = ES ia . 



(18) 
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In addition to the pole part, the S-wave amplitudes 
receive also a contribution coming from the commutator 
in fT^ 39] This is usually referred to as the "soft-pion" 
term: 

A a(sof t) = Z V2 {BMiHeffmy (1Q) 

Unlike the factorization part, the pole and soft-pion 
terms involve matrix elements which are not directly re- 
lated to experiments and have to be estimated theoret- 
ically. In the next section we present our calculation of 
these matrix elements in the ILM. 



IV. ILM CALCULATION 

In this section we present our calculation of the P-wave 
and S-wave amplitudes, within the ILM. 

A. P-wave amplitudes 

In order to determine the P-wave amplitudes in the 
ILM model, we need to evaluate the non-perturbative 
inputs h^, defined in (|14|l . 

In a field-theoretic framework, these matrix elements 
can be extracted from appropriate ratios of Euclidean 
three- and two- point functions. Let us consider the 
three-point correlator: 

Gi' B (r) = Jd^J d 3 y(0|T( J£,(x,2r) 

We//(y,r)^(0,0))|0), (20) 

where r = i t, a is a spinor index and J B (x) , J B , (x) 
are interpolating operators which excite states with the 
quantum numbers of the B and B' baryons. (For 
example, for the proton and S + hyperon we used 
Jp{x) = e abc (uj (x) C 75 d b (x)) u% (x) and Jg+(a;) = 
tabc (s a r (x)Cj 5 u b (x)) v%{x).) 

It is straightforward to show that, in the limit of large 
Euclidean time separation, the correlator l|20(l relates di- 
rectly to the matrix element h B , B : 

lim G? B (r) =2h + B , B A / B A B e-^' +M ^\ 

T — >00 

(21) 

where Ab' and Ab are the couplings of the interpolating 
fields Jb' and Jb' to the B 1 and B states, defined as 

(0\J B (x)\B) = A B u B (p)e^ x . (22) 

In the SU / (3) symmetric limit we are considering we have 
A B > = Ab = A and Mb> = Mb = M. Hence, in this ap- 
proximation, it is possible to extract the matrix clement 
h B , B by taking the ratio of the three-point function 120|) 
with -say- the proton two-point function: 

f<B' Bl \ 



where 

G 2 (r) = y>x(0| T[Jg(x,r)J£(0,0)] |0) 
T ^°° 2A 2 e- Mr . 

(24) 

Non-perturbative calculations of QCD correlation 
functions can be performed by exploiting the analogy be- 
tween the Euclidean generating functional and the par- 
tition function of a statistical ensemble. In lattice QCD, 
one usually carries-out analytically the integral over the 
fermionic fields, and then computes numerically Monte 
Carlo averages of the resulting Wick contractions over a 
statistical ensemble of gauge configurations. In the ILM, 
we replace the space of all gauge configurations with an 
ensemble of instantons and anti-instantons Q. Like in 
lattice QCD, in each configuration the quark propagator 
is obtained by inverting the Dirac operator. Unlike in lat- 
tice QCD, in the ILM there is no need of regularization, 
so all calculations are performed in the continuum. This 
prescription is equivalent to computing the correlation 
functions to all orders in the 't Hooft interaction. 

In this work we have considered the simplest version 
of the model, the Random Instanton Liquid (PJLM), in 
which the density and size of the pseudo-particles are 
kept fixed, while their position in a periodic box and their 
color orientation are generated according to a random 
distribution. 

We have evaluated numerically |3(j the correlation 
functions associated to the matrix elements (p|iJ e ^/|E + ) 
and (A\H e f j\E°) . We have averaged over 52 configura- 
tions of 252 pseudo-particles of size p — 0.33 fm, in a peri- 
odic box of volume (3.6 3 x 5.4) fm 4 . Like in lattice simu- 
lations, we have chosen a rather large current quark mass 
for u and d quarks (75 MeV) , to avoid finite- volume arti- 
facts. In order to check for the dependence of our results 
on the quark masses, we have also performed the same 
calculation using larger quark masses (135 MeV). Finally, 
to enforce flavor symmetry, we have set m s = m u = rrid- 
The 6-dimensional spatial integration in l|20|l has been 
performed by means of an adaptive Monte Carlo method 
(VEGAS). Convergence has been achieved using 1600 in- 
tegration points. The 3-dimensional integral in l|24|) has 
been performed by first carrying out the angular integra- 
tion analytically (exploiting rotational symmetry) and 
then computing the remaining 1-dimensional radial in- 
tegration by a Gauss-quadrature method. 

We have observed that the quark-model relation 
f + /d + ~ 1 holds also in our field-theoretic approach 
gg, with d+ = (0.28 ± 0.05) x 10" 7 GeV, a re- 
sult quite close to the prediction of the SNXD model 
( d+ = 0.35 x 10~ 7 GeV |2jJ ). This calculation 
shows explicitly that gluon and sea degrees of freedom 
contribute very little to these decay amplitudes. 

The results presented so far correspond to sim- 
ulations performed with quark masses of 75 MeV. 
We have found that calculations with heavier quark 
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P-Wave Amplitudes ( x 10 7 ) 



Pole Fact. RILM Experiment 



Ao 


-6.87 


-4.03 


-10.9 ± 1.17 


-15.61 ± 1.4 


0.7 


A 


9.72 


8 


17.71 ± 1.66 


22.40 ± 0.54 


0.8 




20.82 


1.65 


22.4 ±3.55 


26.74 ± 1.32 


0.8 


K 


31.84 





31.84 ±4.81 


41.83 ±0.17 


0.8 




1.75 


-3.26 


-1.52 ±0.30 


-1.44 ±0.17 


1.1 




16.15 


-2 


14.15 ± 2.75 


17.45 ±0.58 


0.8 


—0 


-11.42 


1.01 


-10.42 ± 1.95 


-12.13 ±0.71 


0.9 



TABLE I: Theoretical prediction and experimental results for 
P-wave amplitudes. Following the standard notation, cor- 
responds to Amp(B Q —►£?' + 7r 9 ). The RILM prediction is 
obtained by adding the pole and factorization contribution. 
Wilson coefficients have been evaluated at the hadronic scale 
H = 1/p = 0.6 GeV, using A-^ = 230 MeV. 



masses (135 MeV) lead to very similar results 
( d+ = 0.27 ± 0.04 x 10~ 7 GeV ). We can there- 
fore conclude that the dependence of these amplitudes 
on the quark mass is very weak. 

It is important to ask whether diquark quasi-bound 
states survive within l/2 + baryons, or if they are melted 
by the interaction with the third quark. To answer, we 
have compared matrix elements obtained from the scalar 
and from the pseudo-scalar part of the diquark operator 
in • We have found that such matrix elements are in- 
deed dominated by the scalar operators in the effective 
Hamiltonian. This is a non-trivial result which represents 
clean signature of the existence of scalar diquark struc- 
tures in the hyperons, in the ILM. On the other hand, it 
also implies that pseudo-scalar diquarks are not present 
in such baryons. Finally, since final-state interaction ef- 
fects are presumably small in this channel |25|. we have 
neglected them. 

Our results for the P-wave amplitudes, obtained by col- 
lecting the factorization and the pole contributions, are 
reported in table ^ and compared to experimental data. 
First of all, we observe that the RILM can reproduce the 
overall body of data on P-wave hyperon decays. All the- 
oretical amplitudes lie within approximatively 20% form 
the experimental results. Note that this discrepancy is 
of the order of the systematic error introduced by the as- 
sumption of SU/(3) symmetry. However, taking a closer 
look, we notice that the central values of the theoretical 
predictions consistently undershoot the experimental re- 
sults (except in one case to be discussed below). This is 
hardly surprising, because in the present calculation, we 
have neglected all confining interactions. 

Finally, we observe that the theoretical prediction for 
the amplitude £~ is the only one overshooting the ex- 
perimental data. This is probably a reflection of the fact 
that this is a very delicate channel, where the factoriza- 
tion and pole terms are of the same order of magnitude 
and have opposite sign. 



B. S-wave amplitudes 

S-wave amplitudes, receive contributions from both the 
pole and the soft-pion part of the PCAC amplitudes. 

The pole part involves mixing of l/2 + baryons with 
l/2~ virtual intermediate states (|14|) . As discussed in 
detail in [2^, in a simple quark-diquark model, in or- 
der for the matrix elements to be non- vanishing one 
needs to assume the existence of 0~ diquark structures 
l/2~ octet baryons. On the other hand, the 't Hooft 
interaction is repulsive in the 0~ channel. While the at- 
traction in the + channel triggers the formation of scalar 
diquarks in l/2 + hyperons contributing to P-wave am- 
plitudes, the repulsion in the 0~ channel prevents the 
formation of pseudo-scalar diquarks, which would show 
up in 1/2 ~ hyperons. Hence, in the ILM, the pole contri- 
bution to S-wave amplitudes is expected to be suppressed 
and we shall neglect it. 

On the other hand, we compute explicitly the soft-pion 
term (|19|l . which arises from the commutator in i|12f) . For- 
sake of definiteness, let us consider the (P ir°\H e f 
S-wave transition. The relevant part of the Q_ operator 
in the effective Hamiltonian can be written in a simplified 
notation as: 

(duy Q+ (us)o- ± {du)l_ (us) + + h.c. (25) 

The soft-pion contribution depends on the commutator of 
the effective Hamiltonian with the axial-charge operator. 
Using current-algebra relationships it is possible to show 
that the commutator of H25|) with if gives the same result 
as the commutator of the operator: 

(duy 0+ (us) + + (du)l_ (us)q- + h.c. (26) 

with the I a operator. Due to the repulsion of the 't Hooft 
interaction in the 0~ diquark channel, the instanton con- 
tribution to the matrix elements of the second term in 
(|26|l between l/2 + states is negligible. On the other 
hand, the matrix elements of the first term in Q26[) relate 
to the / + and d + constants, which have been calculated 
to determine the P-wave amplitudes. 

Final-state interaction corrections in this channel are 
rather small but not negligible. We have included them 
following the estimate performed in [24|. 

The RILM predictions for S-wave decay amplitudes are 
presented in table ^ and compared to experimental re- 
sults. As in the case of P-wave transitions, we observe a 
good agreement with experiment, with prediction within 
20% from the data. Again, we observe that the ILM tends 
to undershoot the measured amplitudes, which confirms 
the idea that roughly 20% of the attraction in the + anti- 
triplet channel comes from confining interactions. Note 
that having neglected pole term leads to very reasonable 
results ( except in one channel, where also all other 
contributions vanish). Clearly, there is no need to assume 
pseudo-scalar diquark structures in l/2~ baryons. 
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S-Wave Amplitudes ( x 10 7 ) 



soft fact. RILM RILM(FSI) Experiment 



A 


— 1.71 


0.2 


— 1.51 


— 1.75 ± 0.34 


—2.36 ± 0.03 


0.7 




A 


2.41 


-0.53 


1.88 


2.25 ±0.57 


3.25 ±0.02 


0.7 


S 0.5 

< 




-4.18 


0.23 


-3.96 


-3.55 ±0.64 


-3.25 ±0.02 


1.1 


K 














0.14 ±0.03 






si 


5.91 


-0.62 


5.29 


4.34 ±0.9 


4.27 ±0.01 


1 






-4.83 


0.61 


-4.22 


-4.22 ±0.82 


-4.49 ± 0.02 


0.9 




—0 


3.41 


-0.22 


3.20 


3.20 ±0.58 


3.43 ±0.06 


0.9 






TABLE II: Theoretical prediction and experimental results 
for S-wave amplitudes. Following the standard notation, 
Af corresponds to Amp(B Q + tt 9 ). The RILM pre- 

diction is obtained by adding the soft-pion and factoriza- 
tion contributions. The results in RILM(FSI) include also 
final-state interaction corrections, as estimated in |25| . Wil- 
son coefficients have been evaluated at the hadronic scale 



H = 1/p = 0.6 GeV, using Aj 



230 MeV. 



V. DISCUSSION 

In the previous section we have shown that the in- 
clusion of instanton-induced effects allows to reproduce 
the overall body of data on non-leptonic hyperon decays. 
We recall that both P-wave and S-wave results have been 
obtained by considering only the contribution of the op- 
erator Q_ in the effective Hamiltonian, which drives only 
transitions with violation of isospin 1/2. Hence, we con- 
clude that the 't Hooft interaction does provide a non- 
perturbative dynamical explanation of the A / = 1/2 
rule. 

An important question to ask is whether one can rule- 
out alternative dynamical mechanisms, which are not 
based on quasi-classical interactions. As already stressed, 
the essential dynamical property which is required in or- 
der to produce an enhancement of A / =1/2 transi- 
tions is an attraction in the scalar anti-triplet + channel. 
Clearly, any model for the microscopic dynamics which 
exhibits a sufficiently strong attraction in this channel 
will produce scalar diquarks |41|. It is nevertheless very 
important to clarify the dynamical origin of these struc- 
tures, whose existence seems to be confirmed by a number 
of independent phenomenological studies (for example, in 
connection with exotic spectroscopy, see |32ll33| L 

In the following we suggest a systematic, model- 
independent procedure to answer the question whether 
quasi-classical topological fields do indeed provide the 
dominant non-perturbative interactions driving diquark 
formation and the A / =1/2 rule. The idea is to eval- 
uate the relevant matrix elements on the lattice and to 
compare the behavior under cooling of the decay ampli- 
tudes and of the string tension. 

The cooling algorithm consists of performing statisti- 
cal averages on different ensembles of gauge configura- 
tions which are closer and closer to the extreme of the 
Euclidean action. This way, the contribution of quasi- 



ISCENARIO 1 
ISCENARIO 2 



0.5 

o/o„ 



FIG. 3: Two different hypothetical scenarios for the behav- 
ior of lattice QCD decay amplitudes under cooling. On the 
x axis, a/ao represents the values for string tension obtained 
after different numbers of cooling steps, normalized to the 
QCD string tension (no cooling). On the y-axis, A/Ao repre- 
sents the ratio of a decay amplitude computed after the same 
number of cooling steps, normalized to its value in QCD (no 
cooling). In SCENARIO 1, the decays are driven by quasi- 
classical interactions and the amplitudes change by roughly 
20% under cooling (ILM prediction). In SCENARIO 2, the 
decays are driven by the confining forces and vanish rapidly 
under cooling. 



classical fields is progressively isolated. It is well known 
that, after few cooling steps, all perturbative fluctuations 
as well as the confining interactions are removed from the 
QCD vacuum. On the other hand, the essential proper- 
ties of light hadrons, such as their masses and point-to- 
point correlators, are seen to change very little. This 
implies that light hadrons are predominantly bound by 
quasi-classical non-confining gauge configurations |(| . 

The main shortcoming of the cooling procedure is that 
it leads to results which intrinsically depend on the ar- 
bitrary number of cooling steps. Due to this problem, it 
is very difficult to make systematic, quantitative state- 
ments. On the other hand, the qualitative observation 
that light hadrons still exist in the absence of confine- 
ment and that smooth, topological structures survive 
even when the string tension is drastically suppressed 
are model-independent facts, in QCD. 

We recall that instantons are smooth, topological 
quasi-classical configurations which bind hadrons but do 
not confine. This observation suggests to study the be- 
havior of the A / =1/2 decay amplitudes as a func- 
tion of the string tension, calculated after each cooling 
step (see Fig. |3J). On the basis of our analysis we pre- 
dict that, if instantons are indeed the leading dynamical 
effect, then the amplitudes should decrease by at most 
20%, as the string tension varies from its physical value 
to nearly zero. On the other hand, if instantons do not 
provide the dominant interaction in these processes, then 
the amplitudes should drastically die out, along with the 
string tension. 
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VI. CONCLUSIONS AND OUTLOOK 

In this work, we have studied the instanton contri- 
bution to non-leptonic weak decays of hyperons. We 
have applied the OPE formalism to separate hard-gluon 
corrections to soft non-perturbative effects and we have 
used the Random Instanton Liquid Model to compute 
the relevant low-energy matrix elements. The connection 
between the matrix elements and the decay amplitudes 
has been established considering the contributions arising 
from both the pole and soft-pion terms in the PCAC re- 
lations and from the factorization part of the amplitude. 
Final-state interaction corrections have been applied to 
S-wave transitions, and have been neglected in P-wave 
transitions. 

We have found that the ILM yields to a good de- 
scription of both P-wave and S-wave decays, providing 
a microscopic explanation for the A 7 = 1/2 rule. In 
this model, the strong enhancement of the transitions in 
which the total isospin is changed by 1/2 units is orig- 
inated by the strong attraction due to the 't Hooft in- 
teraction in the quark-quark scalar anti-triplet channel, 
leading to a quark-diquark structure in the hyperons. 
We stress that the calculation presented in this work 
were performed with no parameter fitting. The only phe- 
nomenological quantities introduced by the ILM are the 
instanton average size and density, which have been fixed 
long ago to reproduce global vacuum properties. 



Our results provide a further confirmation of the gen- 
erally accepted picture according to which the internal 
dynamics of light hadrons is dominated by the interac- 
tions responsible for chiral symmetry breaking. Indeed 
in the present calculation, roughly 70% of the ampli- 
tudes comes from instanton-induced interactions (which 
drive the spontaneous breaking of chiral symmetry), 10% 
from hard gluon-exchange corrections, while the remain- 
ing 20% is due to some other interactions, presumably 
related to confinement. Results are seen to depend very 
weakly on the value of the current quark masses chosen. 

Since the present analysis is affected by some model- 
dependence, we cannot in principle rule-out possible al- 
ternative dynamical mechanisms for scalar diquark for- 
mation. However, we have suggested a lattice-based pro- 
cedure which would allow to determine, in a unambiguous 
and model-independent way, if the strong attraction in 
the diquark channel is generated by quasi-classical gauge 
configurations or if it is due to the quantum fluctuations 
associated with the dynamics of color confinement. 
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